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FOURIER  SERIES  OPERATING  PACKAGE 

by 

Milton  L.  Charnow 
Goddard  Space  Flight  Center 


SUMMARY 

This  report  presents  a  computer  program  for  multiplying, 
adding,  differentiating,  integrating,  "barring"  and  scalarly 
multiplying  "literal"  Fourier  series  as  such,  and  for  extracting 
the  coefficients  of  specified  terms. 
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FOURIER  SERIES  OPERATING  PACKAGE 


by 

Milton  L.  Charnow 
Goddard  Space  Flight  Center 


INTRODUCTION 

The  Hansen  Satellite  Theory  as  modified  by  Musen^  involves  various  manipulations 
of  "literal"  Fourier  series,  as  such,  before  final  numerical  evaluation.  To  program  this 
formulation  it  is  necessary  to  represent  cosines,  sines,  and  constants  in  a  special  man¬ 
ner  and  thus  manipulate  the  various  Fourier  series. 


SPECIAL  REPRESENTATION  OF  A  SINE,  COSINE, 

'  OR  CONSTANT 

■i'  It  was  decided  to  represent  any  term,  including  the  constant  term,  of  the  Fourier 

series  we  are  concerned  with  by  two  consecutive  eight-place  floating-point  words.  The 
terms  of  this  series  have  the  general  form 

cos  (iF  +  jE  +  kW  +  lU) 
or 

Ajj  sin  (iF  +  jE  +  kW  +  lU) 

and  the  values  of  F,  E,  W,  and  U  are  not  used  until  final  numerical  evaluation.  The  first 
eight-place  floating-point  word  represents  the  coefficient  An  in  an  entirely  standard  man¬ 
ner.  The  second  eight-place  floating-point  word  represents  the  sine  or  cosine  term  of  up 
to  four  arguments  in  an  artificial  manner. 


^Musen,  P.,  "A  Modified  Hansen's  Theory  as  Applied  to  the  Motion  of  Artificial  Satel¬ 
lites,"  NASA  Technical  Note  D-492,  November  I960. 
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Consider  a  term  such  as  An  cos  (iF  +  jE  +  kW  +  lU).  Until  this  term  is  evaluated  we 
are  concerned  only  with  the  values  of  An,  i,  j,  k,  and  1.  The  coefficient  An  is  represented  # 

by  the  first  of  the  two  eight-digit  floating-point  words.  The  four  coefficients  of  the  argu¬ 
ments  (i,  j,  k,  and  1)  are  each  represented  by  two  digits  of  the  second  eight-place  floating¬ 
point  word.  The  first  two  digits  represent  i  and  are  normalized  to  00;  the  next  three  pairs 
of  digits  represent  j,  k,  and  1,  and  are  normalized  to  50.  Thus,  the  value  of  i  may  vary 
from  0  to  99,  while  j,  k,  and  1  may  vary  from  -49  to  +49.  A  cosine  is  denoted  by  a  plus  C 

sign;  a  sine  by  a  minus  sign,  that  is,  1  cos  (OF  +  OE  +  OW  +  OU)  or  cos  (0)  becomes 
+10000000  +  01,  +00505050  +  08.  Any  constant  term  can  be  represented  as  A^  cos  (0).  A 
few  additional  examples  will  be  helpful  (see  also  Appendix  A): 


Conventional 


Special 


1/4  cos  (IF  +  2E  -  3W  -  2U) 
1/4  sin  (OF  +  IE  -  OW  +  2U) 
1/4 
0 


25000000  +  00, 
25000000  +  00, 
25000000  +  00, 
00000000  +  00, 


+01524748  +  08 
-00515052  +  08 
+00505050  +  08 
+00505050  +  08 


The  convention  of  having  the  first  non-zero  coefficient  of  the  argument  word  positive  was 
adopted.  Cos(-x)  becomes  cos(x)  and  sin(-x)  becomes  -sin(x).  Examples  are: 

(1/2)  cos  (OF  -  3E  +  IW  +  OU)  becomes  (1/2)  cos  (OF  +  3E  -  IW  +  OU) 
and  the  special  representation  is  50000000  +  00,  +  00534950  +  08 

(1/2)  sin  (OF  -  3E  +  IW  +  OU)  becomes  -(1/2)  sin  (OF  +  3E  -  IW  +  OU) 
and  the  special  representation  is  -50000000  +  00  -  00534950  +  08. 


The  first  location  address  of  a  series  contains  the  number  of  terms  of  the  series. 
A  series  of  n  terms  would  be  represented  by  2n  +  1  words,  the  first  of  which  would  be 
the  number  n. 


THE  FOURIER  OPERATING  PACKAGE 


The  following  series  operations  are  performed  by  the  Fourier  Operating  Package: 

•  Multiplication 

•  Addition  and  Subtraction 

•  Differentiation 

•  Integration 

•  Bar  (Special  operation  used  in  the  Hansen  satellite  theory) 

•  Scalar  Multiplication 

•  Coefficient  Extraction 

•  Series  Evaluation 
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Multiplication 

Multiplication  of  two  series,  where  the  terms  are  of  the  general  form  described 
earlier  and  the  values  of  F,  E,  W,  and  U  are  not  used  until  the  final  numerical  evalua¬ 
tion,  is  according  to  the  conventional  trigonometric  identities: 

A  cos  X-B  cos  Y  =  (AB/2)  cos  (X+Y)  +  (AB/2)  cos  (X-Y) 

A  cos  X-B  sin  Y  =  (AB/2)  sin  (X+Y)  -  (AB/2)  sin  (X-Y) 

A  sin  X-B  cos  Y  =  (AB/2)  sin  (X+Y)  +  (AB/2)  sin  (X-Y) 

A  sin  X-B  sin  Y  =  (AB/2)  cos  (X+Y)  +  (AB/2)  cos  (X-Y). 

The  Fourier  Multiplication  routine  is  composed  of  three  major  sections:  the  multiplier, 

the  collapser,  and  the  arranger. 


The  Multiplier 


Two  Fourier  series  such  as 


and 


AiAi  +  A2A2  +  AgAg  +  .  .  .  +  Aj^Ajj 

+  ®2®2  ^  +  •  •  -  +  ®m®m 


(Series  A) 
(Series  B) 


which  are  to  be  multiplied  are  arranged  in  descending  order  of  the  absolute  values  of  the 
coefficients,  that  is, 


Al 


> 


> 


>.  .  .  > 


An 


and 


Bl 


> 


B2 


> 


B3 


>  .  .  .  > 


B 


m 


To  facilitate  further  discussion,  we  shall  denote  any  term  in  the  A  series  AxA|, 

*  * 
any  term  in  the  B  series  as  B  B  ,  and  any  term  in  the  resultant  series  by  C  C„,  where 

y  y  :jt  ♦  z  z 

Aj^,  By  and  are  the  coefficient  words  and  A^,  By,  and  are  the  argument  words. 

The  multiplication  of  the  A  series  by  the  B  series  proceeds  as  follows:  The  first 
term  in  the  A  series  is  multiplied  by  each  term  in  the  B  series,  then  the  second  term  in 
the  A  series  is  multiplied  by  each  terra  in  the  B  series,  and  so  on  until  each  term  in  the 
A  series  has  been  multiplied  by  each  term  in  the  B  series.  For  example^|AxBy|  is  compared 
with  some  numerical  criterion  e.  If  j  AxBy|  >  e,  then  and  +  i)  +  1)  are  gen¬ 
erated  according  to  the  trigonometric  formulas  all  eady  stated.  If  j  AxByj  <  e,  then  the 
value  of  y  is  examined.  If  y  >  1  (i.e..  By  is  any  term  other  than  the  first  term).  Ax  is 
replaced  by  ^  and  A^j^  ^  is  multiplied  by  Bj^.  If  y  =  1,  the  multiplication  of  the 
two  series  is  terminated  since  any  further  A^jj  ^  By  will  also  be  less  than  e.  The  mul¬ 
tiplication  process  continues  until  A^Aj^  has  been  multiplied  by  BjjjBj^  unless  the 
numerical  criterion  or  space  limitations  intervene. 
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The  Collapser 

Every  multiplication  generates  two  terms  of  two  words  each.  The  purpose  of  the 

3|( 

collapser  is  to  combine  like  argument  terms.  Each  argument  term  is  compared  with 

each  other  argument  term  previously  generated  and  stored.  K  C_  equals  any  other  ' 

z 

argument  term,  the  corresponding  coefficient  terms  are  added.  Thus,  there  is  no  dupli¬ 
cation  of  terms, 

The  Arranger 

The  final  step  in  the  multiplication  is  the  arranging  of  the  terms  of  the  series. 

I  Cl  I  is  compared  with  |  C2I,  [C3I,  etc.  If  |C2;  |  >  |  Ci|,  then  Cj  is  replaced  by  and 
Cl  is  replaced  by  C*.  The  process  continues  until  the  terms  are  arranged,  in  descending 
order,  according  to  the  absolute  value  of  the  coefficients. 


Addition  and  Subtraction 

Addition  or  subtraction  of  two  Fourier  series  is  primarily  a  process  of  com¬ 
paring  argument  terms  and  adding  the  coefficients  of  like  terms.  aJ  is  compared  suc¬ 
cessively  with  Bi  through  A2  with  Bf  through  B^,  etc.,  until  has  been  compared 
with  B^.  If  A^  =  By,  the  sum  of  the  coefficients  (Ajj  +  By)  and  the  argument  term  A^  are 
stored,  and  By  and  By  are  replaced  by  zeros  If  A^  does  not  equal  any  By,  both  A^  and 
A*  are  stored.  After  all  terms  in  the  A  series  have  been  compared  with  all  terms  in  the  ^ 

B  series,  the  remaining  B  series  terms  are  stored. 

Subtraction  is  accomplished  in  like  manner  after  changing  the  signs  of  each  coeffi¬ 
cient  term  in  the  B  series.  The  resultant  series  in  each  case  is  processed  through  the 
arranger. 


Differentiation 


Differentiation,  in  this  application,  is  with  respect  to  the  F  variable.  Thus, 
A  sin  (iF  +  jE  +  kW  +  lU)  =  iA  cos  (iF  +  jE  +  kW  +  lU). 

Example: 

sin  (3F  +  2E  -  3W  +  U)  =  +3  cos  (3F  +  2E  -  3W  +  U) 


[10000000  +  01,  -03524751  +  08  becomes  +30000000  +  01,  +03524751  +  08. 

On  completion  of  the  differentiation,  the  resultant  series  is  processed  thru  the 
arranger. 


Integration 

Integration,  in  this  application,  is  with  respect  to  E.  However,  W  is  also  a  func¬ 
tion  of  E.  Thus, 

fA  cos  (iF  +  jE  +  kW  +  lU)dE  =  ■  .  sin  (iF  +  jE  +  kW  +  lU). 

J  Cj  ]  +  C2  K 

Example : 

With  C|  =  1  and  C2  =  1, 

j[30000000  +  01,  +  01525250  +  08]dE  becomes  75000000  +  00,  -01525250  +  08. 

The  resultant  integrated  series  is  also  processed  through  the  arranger. 


Bar  Operation 


The  bar  operation  is  a  special  function  in  the  Hansen  Satellite  Theory.  It  con¬ 
sists  of  adding  the  coefficient  of  the  F  argument  to  the  coefficient  of  the  E  argument  and 
substituting  zero  for  the  F  coefficient.  Thus, 

A  cos  (iF  +  jE  +  kW  +  lU)  after  barring  becomes  A  cos  (OF  +  (i  +  j)  E  +  kW  +  lU). 

Example : 

50000000  +  00,  +  02534850  +  08  after  barring  becomes  50000000  +  OC,  +  00554850  +  08. 


Scalar  Multiplication 

Scalar  multiplication  is  the  multiplication  of  the  coefficient  A^^i  of  each  term  by 
a  constant. 


Coefficient  Extraction 


In  the  Hansen  Satellite  Theory  it  is  sometimes  necessary  to  use  the  coefficient 
of  some  term  of  a  Fourier  series  such  as  a  sine  IF  term,  cosine  2E  term,  or  the  con¬ 
stant  term  of  a  series.  Let  us  assume  it  is  necessary  to  use  the  constant  term.  If  there 
is  a  constant  term  in  that  series,  it  will  be  the  multiplier  of  the  cos  (0),  or  in  special 
form,  that  word  which  multiplies  00505050  +  08.  We  successively  compare  each 
argument  term  in  the  series  with  cos  (0)  and  extract  that  which  multiplies  the  argu¬ 
ment  term  cos  (0).  If  no  argument  term  of  the  series  is  cos  (0),  a  normalized  zero 
(10000000  +  00  +  00000000  +  00  +  00505050  +  08),  is  printed. 


Series  Evaluation 

The  numerical  values  F,  E,  W,  and  U  are  only  employed  in  the  Series  Evalua¬ 
tion  Routine. 

To  evaluate  a  Fourier  series,  the  numerical  values  of  i,  j,  k,  and  1  are  multi¬ 
plied  by  the  numerical  values  of  F,  E,  W,  and  U,  and  the  sum  iF  +  jE  +  kW  +  lU  is  de¬ 
termined.  The  sine  or  cosine  of  iF  +  jE  +  hW  +  lU  is  multiplied  by  the  coefficient  A  and 
the  terms  are  added. 


CONCLUDING  REMARKS 

Appendix  A  is  the  special  representation  of  a  nine  term  series.  Appendix  B  presents 
flow  charts  and  Appendix  C  a  listing  of  instructions  for  the  program  packages  for  multi¬ 
plication  (including  the  collapser  and  arranger),  addition  or  subtraction,  differentiation, 
integration,  bar  operation,  scalar  multiplication,  coefficient  extraction,  and  the  evalua¬ 
tion  of  the  final  series.  Because  this  program  was  written  in  Mystic  Code  for  the  IBM 
709,  an  explanation  of  Mystic  Code  is  given  in  Appendix  D. 

The  Fourier  operating  package  can  be  used  with  any  theory  that  involves  representa¬ 
tions  of  functions  by  Fourier  series.  It  can  also  be  modified  to  operate  with  polynomials 
of  the  form 
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Appendix  A 

Nine  Term  Series  in  Standard  and  Special  Form 

Q  The  following  is  a  nine  term  series  presented  in  standard  and  in  special  form.  Note 

that  in  the  special  form,  the  first  location  (address)  contains  the  number  of  terms  in  the 
series. 


Series  in  standard  representation 

.29467121  cos  (0) 

+.00010496334  cos  (OF  +  OE  +  2W  +  OU) 
+.00005252596  sin  (OF  +  lE  -  IW  +  OU) 
+.000019845618  cos  (OF  +  IE  -  2W  +  OU) 
+.0000066329604  cos  (OF  +  IE  +  2W  +  OU) 
+.0000020107054  cos  (OF  +  2E  -  2W  +  OU) 
-.00000036004597  sin  (OF  +  OE  +  IW  +  OU) 
-.000000055052357  sin  (OF  +  OE  +  3W  +  OU) 
-.000000031090653  sin  (OF  +  IE  +  IW  +  OU) 


Series  in  special  representation 

90000000  +  01 

+29467127+00,  +00505050+08 
+10496334-03,  +00505250+08 
+52525962-04, -00514950+08 
+19845618-04,  +00514850+08 
+66329604-05,  +00515250+08 
-20107054-05,  +00524850+06 
-36004597-06,  -00505150+08 
-55052357-07,  -00505350+08 
-31090653-07, -00515150+08 
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Appendix  B 

Flow  Charts 

The  following  are  the  flow  charts  for  the  op¬ 
erating  packages  for  multiplication  (including 
the  collapser  and  the  arranger),  addition  or 
subtraction,  bar  operation,  differentiation,  in¬ 
tegration,  and  series  evaluation. 
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Flow  Chart  for  Multiplication  (Continued) 
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Flow  Chart  for  Multiplication  (Continued) 


Coliapser 
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Fl^w  Chart  for  Bar  Operation 


Flow  Chart  for  Differentiation 


Initiate  2N 
Initiate  a 


mm 


u 


/  Determine  \ 
A,  =  A^/kg 
^as  A,  sin  or  cos? 


Appendix  C 


Listing  of  Instructions 

The  following  is  a  listing  of  instructions  for 
multiplication  (including  the  collapser  and  the 
arranger),  addition  or  subtraction,  differentia¬ 
tion,  integration,  bar  operation,  scalar  multi¬ 
plication,  and  coefficient  extraction. 

The  square  root  instructions  are  for  stand¬ 
ard  number  representations  and  therefore  have 
not  been  discussed  in  the  report. 
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Appendix  D 


Mystic  Code 

The  Fourier  Operating  Package  described  in 
this  report  was  written  for  the  IBM  709  in  Mys¬ 
tic  Code.  The  following  is  an  explanation  of  the 
Mystic  Code. 
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